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We consider the problem of evaluating the entanglement of non-Gaussian mixed states generated 
by photon subtraction from entangled squeezed states. The entanglement measures we use are 
the negativity and the logarithmic negativity. These measures possess the unusual property of 
being computable with linear àlgebra packages even for high-dimensional quantum systems. We 
numerically evaluate these measures for the non-Gaussian mixed states which are generated by 
photon subtraction with on/off photon detectors. The results are compared with the behavior of 
certain operational measures, namely the teleportation fidelity and the mutual information in the 
dense coding scheme. It is found that all of these results are mutually consistent, in the sense that 
whenever the enhancement is seen in terms of the operational measures, the negativity and the 
logarithmic negativity are also enhanced. 



PACS numbers: 03.67.Mn, 03.67.Hk, 42.50.Dv 



I. INTRODUCTION 

Continuous variable (CV) quantum optical systems are 
well-established tools for both theoretical and experi- 
mental investigations of quantum information process- 
ing (QIP) QjQ- The essential resource, entanglement, 
can be realized with a Gaussian two-mode squeezed vac- 
uum state. This state is relatively easy to work with 
theoretically and is also commonly produced in the lab- 
oratory. It has been successfully applied to implement 
various important protocols, such as quantum teleporta- 
tion d, 0, IS @ i quantum dense coding 0, H, @ and en- 
tanglement swapping 0, . These successes are based 
on well-developed techniques of optical Gaussian oper- 
ations. These consist of beam splitting, phase shifting, 
squeezing, displacement and homodyne detection. 

However, recent theoretical investigations have shown 
some limits on such operations. A prime example is the 
no-go theorem relating to the distillation of entanglement 
shared by distant parties using only Gaussian local oper- 
ations and classical communication (LOCC) p^ . lÏ3lll^| . 
To go beyond this limit, one should use higher order 
nonlinear processes, such as the cubic-phase gate |Ï5j or 
Kerr nonlinearity |16(. It is, however, difficult to imple- 
ment these nonlinear processes with presently available 
materials, which do not have sufficiently high nonlinear- 
ity and suffer from losses. A more practical alternative 
method has been developed, which uses nonlinear pro- 
cesses induced by photon counting on tapped-off beams 
from squeezed states 0, 0, When ideal photon 
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number resolving detectors are used, the input Gaussian 
state can be transformed into a non-Gaussian pure state 
with higher entanglement. In practice, however, such de- 
tectors are not yet suitable for practical use. The most 
reliable type of photodetector available at present is the 
on/off type photon detector based on avalanche photodi- 
odes. This device can only distinguish the vacuum ('off') 
state from non-vacuum ('on') states. The latter events 
result in a non-Gaussian mixed state. Evaluating the 
entanglement of such a state is far from trivial. 

Previously, the effect of entanglement has been theo- 
retically analysed based on the figures of mèrit of concrete 
protocols, such as the fidelity of teleportation [2CÍ. the de- 
gree of violation of Bcll-type inequalities yïl [22, |23] and 
the mutual information of dense coding |24| . In fact, it 
was shown that the performance of every protocol was 
improved, implying that the entanglement of the non- 
Gaussian mixed state must be enhanced. However, these 
indirect evaluations were dependent upon some external 
parameters, which depend upon specific operations dif- 
fering from protocol to another, such as the type of input 
state for the teleportation fidelity, the choice of the mea- 
surement basis for the Bell-type inequality violation, and 
the signal power of modulation for the dense coding. 

Quantifying the entanglement of the photon- 
subtracted squeezed state, in a way which is independent 
of particular external parameters, is our main concern in 
this paper. To our knowledge, this problem has not been 
previously addressed. For a state p, an entanglement 
measure E(p) should satisfy the following criteria [2o| : 
(i) E is the non-negative functional, (ii) E vanishes if the 
state p is separable and (iii) E should not increase on 
average under LOCC. In general, a quantity satisfying 
these criteria is known as an entanglement monotone. 
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Many such quantities have been proposed, such as the 
entanglement of formation [26L |27| . the entanglement 
cost > the distillable entanglement |2^|, and the 
relative entropy of entanglement |29|. It is, however, 
not easy to calculate these measures for genèric mixed 
states. Recently, however, two entanglement measures 
which are much more amenable to evaluation have been 
proposed. These are the negativity and the logarithmic 
negativit y 1301 . These measures are based on the Peres 
criterion |3Ïj. That is, they are defined in terms of the 
eigenvalues of the partially-transposed density operator. 
The most distinctive feature of these entanglement mea- 
sures is that they are easily computable numerically with 
linear àlgebra packages. Furthermore, the logarithmic 
negativity is an additive functional, and it is an upper 
bound on the distillable entanglement En [3(ï| . 

In addition to being an entanglement monotone, it was 
believed that an entanglement measure E should also be 
(downward) convex, i.e. it should be non-increasing on 
average under the loss of classical information by mixing. 
Convexity is necessary for an entanglement measure to 
be bounded from above by the entanglement of formation 
Ef H2 , although the logarithm function is concave (up- 
ward convex) . Recently, however, it has been shown by 
Plenio et al [3^, |^ that convexity is not directly related 
to the physical process of discarding quantum informa- 
tion, i.e. discarding subsystems such as local ancillas 
during LOCC operations and that the logarithmic neg- 
ativity is indeed a full entanglement monotone |34|. In 
view of these considerations, it is interesting to evalu- 
ate the entanglement of photon-subtracted mixed state 
in terms of the (logarithmic) negativities, and to compare 
them the figures of mèrit associated with the aforemen- 
tioned protocols. Of particular interest to us is how the 
entanglement, and these figures of mèrit, are enhanced 
by the photon-subtraction procedure. 

This paper is organized as follows. In Sec. [HJ we briefly 
summarize the measurement-induced non-Gaussian op- 
eration on the two-mode squeezed vacuum state and its 
mathematical description. In Sec. IIIII the negativity 
and the logarithmic negativity are briefly reviewed, and 
their monotonicity is discussed. In Sec. II VI our mi- 
merical methods for calculating the negativities of the 
photon-subtracted mixed state are presented. In Sec. [V] 
we review the previous analyses of the entanglement of 
such states using protocol-specific figures of mèrit. We 
compare these results with those relating to the nega- 
tivity/logarithmic negativity. The final section IVÏ1 is dc- 
voted to discussion and conclusion. 
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FIG. 1: Measurement-induced non-Gaussian operation on 
the two-mode squeezed vacuum state. BS, PS are beam split- 
ter, photon detector, respectively. 



single-mode squeezed vacuum states, 

\r) k = S k (r)\0) k , (1) 
where Sk(r) is the squeezing operator at path k, 

§ k (r)= eX p[-^-à 2 k )], (2) 

and r is the squeezing parameter. These are combined 
via a balanced beam splitter to generate the two-mode 
squeezed vacuum state, 
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where 
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(3) 
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is the beam splitter operator, and the parameter 9 is 
related to the transmittance T as 



tan# 
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and 9 — tt/A corresponds to the balanced beam splitter. 
The operator Sï, (r) is the two-mode squeezing operator, 



II. MEASUREMENT-INDUCED 
NON-GAUSSIAN OPERATION 

The schematic of the measurement-induced non- 
Gaussian operation on the two-mode squeezed state is 
shown in Fig. ^ The primary sources are two identical, 



^lf( r ) = ex P ~ r ( à Wi - a kài) 



(6) 



Introducing A = tanhr, the Schmidt coefficients aregiven 
by 



yl — A 2 A™. 



(7) 
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The beam at path C (D) is then tapped off from path A 
(B) by a beam splitter of transmittance T. The resulting 
four-mode state just after the second beam splitters is 

IV 1 ) ABCD = CD 
n 

= E a ™ E £mt,nj\n-i)A\n- j)v\ï)c\Ú)v>, 

(8) 

where 

6 lfe = (-i) fe ^Q (Vr)"^(Vïï) fc , (9) 

and (?) is the binomial coefficient and R = 1 — T is the 
refiectance. 



A. Photon number resolving detector case 

When k photons are detected in the beam at path C 
and l are detected in path D by ideal photon number 
resolving detectors, the conditional state is given by 

IV'ngUb OC d(Z|c(*#)aBCD 



E Ot n £, nk (,nl\n - k) \\n - í) B , 
n-max{fc,l} 



(10) 



where |V'ng)ab is still a pure state. In the case of k 
1 = 1, 
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where 
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(12) 



is the probability of detecting one photon in each arm. 
The state 1)110 is not Gaussian any more. 



B. On/off type detector case 

An ideal on/off type detector is described by a positive 
operator-valued measure (POVM) with elements 



rï(° ff ) = 

n(°") - |l)(l| + |2)(2|+··· = ï-|0)(0|. 



(13) 



The two-mode squeezed state is transformed into a mixed 
non-Gaussian state 



Tr C D 



PNG 



' (ABCD) < 



( ft 



(on) 
C 



n 



(on) 



*tí>(AB)<*«l. 



where 

l$ij>AB 



VA\n-j) 



(14) 



(15) 



n-max 



and Vdot is the probability of detecting at least one pho- 
ton in each of the paths C and D, 
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The mean photon number of the state /5ng i s expressed 



as 
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In Fig. |2 the mean photon number Ang is shown as 
a function of A, with that of the two-mode squeezed vac- 
uum state, Asq = 2A 2 /(1 — A 2 ), for comparison. The 
transmittances of the tapping beam splitters are chosen 
as T = 0.9. 

The increase of the mean photon number for the 
photon-subtracted squeezed state is due to the fact that 
the generation process is based on the event selection 
of the components of higher numbers of photons by ex- 
cluding the original vacuum component of the input two- 
mode squeezed state. 



III. COMPUTABLE ENTANGLEMENT 
MEASURES 

In this section, we briefly review the negativity and 
the logarithmic negativity as computable entanglement 
measures that possess the properties of an entanglement 
monotone |30l |: (i) The entanglement measure E is a 
non-negative functional, E(p) > 0, (ii) if p is separable, 
E(p) — 0, and (iii) E(p) does not increase on average 
under the LOCC. 

The negativity of a bipartite mixed state p, denoted 
by N{p) , is defined as the absolute value of the sum of 
the negative eigenvalues of p PT , the partial transpose of 
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FIG. 2: Mean photon number of the photon-subtracted 
mixed state (thick sòlid line) and original two-mode squeezed 
vacuum state (dotted line). 



p with respect to either subsystem. We may write this 
as 



1 



, (18) 



where 1 1 • 1 1 denotes the trace-norm. It is quite easy to 
prové that Aí>0. For a separable state, p = J^íPíPaí <8> 
Pbí, the partial transpose with respect to either subsys- 
tem, say B, is given by p = ^ PíPm <£> p Bi - This is also a 
state, and it therefore has zero negativity. Furthermore, 
under LOCC, M does not increase on average [33,113 ■ ^ 
is also convex, i.e. 

$>AA(À) < Aí(p), (19) 

i 

where p — YliViPi- The logarithmic negativity is defined 
E N {p) = log 2 (l + 2N{p)) = log 2 (||p FT ||). (20) 



as 



In addition to the properties (i), (ii), and (iii), this quan- 
tity is additive because | \p PT <g> a PT \ \ = \ \p PT \ \ ■ \ \Ò PT \ |, 
which is also a desired property of a good entanglement 
measure. 

For a pure entangled state 



Ix)ab = ^2c n \n) A \n) B , 



(21) 



the negativity and the logarithmic negativity can be cal- 
culated analytically, where, without loss of generality, we 
take the Schmidt coefScients c„ to be non-negative. We 
use the fact that 



{|X)(AB)(X|} 



PT 



(22) 



EjïQx)) = 21og 2 (^Ç c «) 



(23) 
(24) 



For example, one can calculate those of the squeezed vac- 
uum state © as 



AA(|r< 2 ))) 



A 



1-A 1 



(25) 



EjsQrW)) = log 2 (l + A)-log 2 (l-A), (26) 

where we have used the Schmidt coefficients described in 
Eq. ©. 

As another example, we here consider the ideal limit 
of photon-subtracted squeezed state (T — > 1). In this 
limit, the photon-subtracted squeezed state with on/off 
detector /5ng is exactly identical to the pure state case 
IV'ngKV'ngIj a lthough the detection probability P^} t ap- 
proaches zero. From Eq. l(TT)l. 
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and thus we have 
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(1 + A 2 )(1-A)' 

log. 
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(27) 

(28) 
(29) 



As a consequence, we find that the logarithmic negativity 
of the photon-subtracted squeezed state is always better 
than that of the original squeezed vacuum for nonzero 
squeezing, i.e. 



E M {\r^)) < limE^l)), 



(30) 



where the equal sign is vàlid for A = 0. 

In realistic situations, on the other hand, the photon 
subtraction would be done by on/off detectors and the 
tapping beam splitters with T < 1. Then the generated 
states are inevitably reduced to mixed states and it is 
typically not possible to obtain analytical expressions for 
their negativity and logarithmic negativity. However, as 
shown in the next section, one can often compute them 
numerically using only linear àlgebra packages. 
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IV. NUMERICAL EVALUATION OF 
NEGATIVITIES 

In this section, we explain the procedure for numerical 
evaluation of the negativities of the photon-subtracted 
squeezed state with the on/ofï detectors. First, we ex- 
pand the output non-Gaussian state (|14J) in the Fock ba- 
sis, 

PNG = E P mi m 2 n 1 n 2 \ 1 nn.l) K{m 2 \ ® B ^2 | , 

m\ ,m2 ,ri2 

(31) 

where 

Pmim2nin2 

= b^iIa^iIpngI^aI^b 

^ oo 

I det , 
«J = l 

(32) 

which means that the density matrix elements are zero 
unless mi — ni = m 2 — n 2 ■ 

The partial transpose of this state, with respect to 
mode B, is 



~PT 

Png 



Pm 1 m 2 n 2 n 1 |mi) A (TO 2 | ® |ni) B (n 2 |, 
m\ ,ï?i2 ,Tii ,ri2 

(33) 

where the elements are zero unless m\ + ni = m 2 + n 2 = 
K. The parameter K is the total photon mimber of both 
beams at paths A and B, and the partially-transposed 
density operator is block diagonal in Fock state basis, 
where the blocks correspond to K = 0, 1, 2, . . .: 



K=0 



(34) 



Here, p^q(K) is the K-th submatrix which is a (K + 1) x 
(K + 1) real matrix (Fig. 0). 

The eigenvalues are obtained by numerically diagonal- 
izing the partially transposed density matrix, 



Q = Ü T p^Ü, 



(35) 



which can separately be done for every submatrix one by 
one, 

K 

= Y,"*\K-l) A {K-l\®\l)A% (36) 

1=0 

where lú? is the Z-th eigenvalue of _ftT-th submatrix, and 

oo 

Ü = Ü K . (37) 

K=0 
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K„ 




FIG. 3: Partially-transposed density matrix of the photon- 



subtracted mixed state with the cut-off K n 
diagonal in Fock state basis. 



which is block 



In the numerical calculation, we introduce a cut-off -Km ax , 
which must be large enough compared with the mean 
photon number. Then, we sum up all the negative eigen- 
values, which gives the negativity, 



1 



NG ) 



A 



K=0 uf <0 



K\ 



(38) 



where 



(39) 



K=0 



is the trace of the truncated density matrix. With this 
result, we can obtain the logarithmic negativity straight- 
forwardly using Eq. 1)20(1 . 

Throughout this paper, we set T = 0.9 unless other- 
wise indicated. In our computations, we chose a cut-off 
of K max — 50. The trace ()39|l could be used as a mea- 
sure of the validi ty of our chosen cut-off. For A = 0.78 
and 0.88, which corresponds to Nng — 7.71 and 14.7, 
A/f max=5 o ~ 1.000 and 0.995, respectively. For the lat- 
ter case, indeed, it still misses 0.5% of the support, but 
this precision sufhces for our analysis. In this connection, 
with -Kmax = 50, it takes 3.5 days to calculate the neg- 
ativity for a certain A, by means of a Mathematica pro- 
gram on Pentium 4 3.2E GHz PC (including the genera- 
tion of matrix components). It is necessary to reach a sat- 
isfactory compromise between the precision of the numer- 
ical calculation and the calculation cost. The numerical 
calculation becomes progressively time-consuming with 
increasing if max . 

Figures 0] and |5] show the numerical vàlues of the 
negativity and the logarithmic negativity, respectively. 
The thick dashed line represents the photon-subtracted 
squeezed state of the pure state case (|ll|l , while the thick 
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FIG. 4: Negativity of the photon-subtracted squeezed states, 
in pure state case (thick dashed line) and mixed state case 
(thick sòlid line). The dotted line corresponds to the input 
squeezed vacuum. 



FIG. 5: The same plot as Fig. 2]of the logarithmic negativity. 
The vàlues of each intersection are A[ N ~ 0.897 and A£? N ~ 
0.772. 



sòlid line corresponds to the mixed one Q14J1. The dot- 
ted line is for the input two-mode squeezed vacuum state 
(|26[l . As can be seen in the figures, the photon sub- 
tracted non-Gaussian states have a larger amount of en- 
tanglement than the input two-mode squeezed vacuum 
state in a practical squeezing range, A < A^ N = 0.897 
for the pure state case, and A < A^ N = 0.772 for the 
mixed one, which correspond to 8.9 dB and 7.1 dB ideal 
squeezing. Exceeding these Aln's, however, the mèrit of 
non-Gaussian operation disappears. This is in contrast 
to the case of T — > 1, where the curve for the photon- 
subtracted squeezed state is always upper than that of 
the squeezed vacuum state, as mentioned in Eq. (I3U|) . 
For A < 0.2, that is, Nng íS 0.28, the difference between 
the cases of the photon number resolving detector and the 
on / off type detector is almost negligible compared with 
the difference between them in the case of ideal two-mode 
squeezed state. This means that the photon subtraction 
by a beam splitter of transmittance T — 0.9 and the 
on/off detectors emulate the single photon subtraction 
well. 



V. OPERATIONAL ENTANGLEMENT 
MEASURES 

In this section, we review the previous results on the 
operational cntanglcmcnt mcasures, and compare them 
with the above results of logarithmic negativity. 



A. CV teleportation fidelity 



(Bell measurement) 




Entangled 
state 



FIG. 6: CV teleportation scheme. BS and HD mean the 
beam splitter and homodyne detection, respectively. 



a classical channel. The schematic is shown in Fig. [ü] 
First, the sender Alice and the receiver Bob share the 

( E'l 

entangled state p AB . Alice then performs the Bell mea- 
surement, that is, a projective measurement in the max- 
imally entangled basis 



|II(a;,p))fcj = 



1 



dyé 



y)k\y)i 



(40) 



upon the unknown state p^ 1 and her fragment of en- 
tangled state (A). She obtains a couple of measurement 
results (x,p), and Bob's part of the initial state (B) is 
correspondingly transformed into 



The first operational measure is the CV teleporta- 
tion fidelity |20j . In quantum teleportation, an unknown 
quantum state can be transferred using entanglement and 



à B {x,p) 



.(in) 

Pt 



.(E) 

Pab 



|n(x,p)) T A 



(41) 
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where 
P (tlp) (x,p) 



Tr B [ TA (n(x lP )\ [fè n) ®p { Aú\ |n(z,p)) T A 

. (42) 

is the Bell measurement probability distribution. Finally, 
Bob corrects his state with the unitary transformation, 
that is, the displacement operation D(pt) — exp[aàT — 
a*à] according to Alice's measurement result, 



x + ip 



ó B (x,p)Dl, 



x + ip 



V2 J"^-"'-» V V2 
The average fidelity between the input and output is 

F = i <l.r I dpP {tlp) (x,p)Ti 



(43) 



p^p {out \x,p) 



(44) 

which depends upon the input state p^ . Let us consider 
that the input state is the coherent state, 



.(in) 



(45) 



which is one of the unaffected state. With the two-mode 
squeezed vacuum state (|3J), the average fidelity is 



^sq(A) 



1 + A 



(46) 



and, with the non-Gaussian pure state (|l·l|l and mixed 
one (THE|) . 



^ng(A) 



(1-A 2 )A 2 T 2 (l-AT+^) 



2P (1) n 



^ng(A) 
respectively, where 



AT) 3 

(A) - ^io(A) - ^oi(A) + ^oo(A) 



1- A 2 



2Pdet 1 - AT — A 2 7í7j - ^(7» + 7i) 



(47) 
(48) 



(49) 



and 71 = R, 70 = 0. See appendices IXI and IBI for morc 
detailed derivations. 



In the limit as T — * 1, .Fng(A) = ^ng(A) and 
F SQ (X) < YimF^(X), 



(50) 



where the equal sign is vàlid when A = and 1. For 
T < 1, however, the average fidelities l|47(l and (|48J) do 
not attain unity even as A — > 1. Figurc \7\ shows the av- 
erage fidelities for T = 0.9. We can see that the photon- 
subtracted squeezed states, in both the pure and mixed 
state cases, are superior to the original squeezed vacuum 
state for A < AÇ = 0.815 and A < A^í = 0.708, re- 
spectively. In the range of AÇ < A, on the other hand, 
the squeezed vacuum state shows the best performance, 
which is quantitatively consistent with the logarithmic 
negativity result. As seen by comparing the results of 
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X 

FIG. 7: Average fidelity of CV teleportation of coherent state, 
with the photon-subtracted pure state (thick dashed line), 
mixed one (thick sòlid line), and original squeezed vacuum 
state (dotted line). The vàlues of each intersection are AÇ ~ 
0.815 and \^ ~ 0.708. 
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FIG. 8: CV dense coding scheme scheme, with power of signal 
modulation a s . BS, M, and HD mean beam splitter, mirror, 
and homodyne detection, respectively. 



Fig. and [7| whenever the fidelity is improved by the 
photon subtraction, so is the logarithmic negativity, i.e. 
At < Aln- We cannot, however, exclude the possibil- 
ity that this is not true for the other input states to 
be teleported. To clarify the exact relation between the 
improvements of the logarithmic negativity and the tele- 
portation fidelity, one has to optimize every component 
of teleportation protocol, such as measurements, etc, for 
the photon-subtracted entangled resource over all pos- 
sible input states. These might be a highly non-trivial 
task. 



B. Mutual information in CV dense coding scheme 

The second operational measure relates to the CV 
dense coding scheme [24j. The schematic is illustrated 
in Fig. [SI In this scheme, Alice and Bob share the entan- 
gled state p A ^ initially. Then, Alice encodes a classical 
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message by a modulation on the beam path A, 

ü A (x s ,p s )=e-W 2 >^D(^±S^j 

(51) 

where x s and p s are related to the power of signal 
modulation a s . For the sake of simplicity, we con- 
sider the quaternary phase-shift keying (QPSK), with 
equal likelihood: a o = (x s = \/2/3,p s = v2/3), 
ara = (V2/3, -V2/9), 010 = (-V2/3, a/2/3), au = 
(-a/2/3, -a/2/3) for some real /3, and P(a fc; ) = 1/4. Bob 
decodes the signals by the Bell measurement Q40p. and 
the decision rule is as follows: {600 = ( x > 0,p > 
Q),6 01 = ( x > 0,p < 0),òi = (x < 0,p > 0),òn = 
(x < 0,p < 0)}. The channel matrix of CV dense cod- 
ing p( ch ) (b mn \aki) > whose elements are the conditional 
probabilities, is calculated with the homodyne probabil- 
ity distribution, 

pW(x,p\x s ,p s ) 

= AB (n(x, p) I Í/a («s , Ps)p AB Ü{( x s i Pa)! 1 ^) P)> AB • 

(52) 

For example, p( ch ) (òoo|aoo) element of the channel ma- 
trix is 



and 



p(ch)( 6( 



ool a ooJ 



dx / dpP 
Jo 



(HD) 



and the other components are calculated similarly. With 
this 4x4 channel matrix, the mutual information is cal- 
culated as 



\a,ki 



J(A;B) = P(aki)P {ch) (b r , 

P {ch) (b mn \a kl ) 



k.L: 



X log 2 



Y Jk >MP{ak'v)P^Hb mn \a k , l ') 



[bit]. 
(54) 



With the two-mode squeezed vacuum state ©, 
/sq(A;B) 



1 + erf i 



1 - erf I 



1 + A 

1 - A 

ï—x 1 



log 2 
log 2 



1 + erf | 



1 - erf 



1 + A 
1 - A 



(3 



1 + A 
1 - A 







where 



erf(ir) = —= \ dte 
\A Ja 



(55) 



(56) 



is the error function. With the non-Gaussian pure state 



4g(A;B) 



-(Jll0g 2 /l +/ 2 l0g 2 /2 

+/ 3 log 2 / 3 +/ 4 log 2 / 4 ), (57) 



h = T>u- 



h = h = 

h = vj- 



1 + erf I 



1 - erf | 



1 - erf I 



1 + AT 



1 - AT 



M/3 



,(58) 



1 + AT 
1 - AT 



/z/3 



(59) 



1 + AT 
1-AT 



/t/3 



,(60) 



where 



(1 - A 2 )A 2 T 2 



P«(1-AT) 2 (1-A 2 T 2 ) 



P 



AT 



2 a2 



2AT a 



1 + AT/ dp? 1 + AT d/i 



,(61) 



containing the differential operation with respect to the 
auxiliary parameter /i. With the non-Gaussian mixed 
state {Q}, 

Tng (A;B) = - (li log 2 li + T 2 log 2 T 2 

+T 3 log 2 T 3 + 14 log 2 J 4 ), (62) 

1 

1 X = ^(-ír+^Il + erf^^)] 2 , (63) 

i,j=0 

1 

£ (-1)^'^ [1 - erf(í^/3) 2 ] , (64) 



where 



T 2 =T: 



1 



and 



li 



C — 



£ (-1)^'^ [1 - erf(^/3)] 2 , (65) 



1-A 2 



4Pdctl-A 2 (T + 7i )(T + 7 J )· 

/ l_ A 2 (r + 7i ) (r + 7 .) 

V (1 — AT) 2 — A 2 7,7j ' 



(66) 
(67) 



where 71 = R and 70 = as before. See appendices IXI 
and IBI for more detailed derivations. 

In the limit as T -> 1, Z NG (A;B) = (A;B) which 
has no intersection with 7sq(A;B) unless A = 0, therefore 



7 SQ (A;B)< üm/«(A;B). 



(68) 



For T < 1, however, above inequality cannot hold for 
larger A. In Figs. lüland lTÜI the mutual information with 
the photon-subtracted squeezed states (T = 0.9) and the 
squeezed vacuum state are shown, at /3 = 1.5 and 0.7, 



9 




P=1.5 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 

FIG. 9: Mutual information in dense coding scheme with the 
photon-subtracted pure state (dashed line), mixed one (sòlid 
line), and squeezed vacuum state (dotted line) at (3 = 1.5. 
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FIG. 11: Intersections of the input squeezed state curve 
and photon-subtracted squeezed state curves. It converges 
to Ad ^ 0.894 for the pure state case, and Aff ~ 0.762 for 
mixed one, as (3 — > 0. 
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X 

FIG. 10: The same plot as Fig. 0at /3 = 0.7. 

respectively. From these results, we can see that the 
range of A for which the non-Gaussian operation brings 
the gain to the mutual information becomes wide as f3 be- 
comes small. When (3 gets smaller, the overlap between 
the probability distributions of different signals increases, 
and the distinction of encoded signals becomes more dif- 
ficult. In this situation, we can see the bona fide cffcct 
of entanglement in CV dense coding scheme. As (3 — > 0, 
the intersections of the curves for the non-Gaussian states 
and the squeezed vacuum state approach Aj^ ~ 0.894 and 
Ap ~ 0.762, respectively (Fig. Illfl . These intersections 
are quite close to those of the logarithmic negativity. This 
is consistent with the logarithmic negativity result in the 
sense that the mutual information indicates the range 
for which non-Gaussian operation enhances the entan- 
glement. 



VI. DISCUSSION AND CONCLUSION 

In this paper, we have studied how to quantify the 
entanglement of the non-Gaussian mixed state which is 



generated by the photon subtraction from the two-mode 
squeezed vacuum state. In order to enhance the entangle- 
ment of a Gaussian state, shared by spatially-separated 
parties, using LOCC, then non-Gaussian operations are 
necessary. When the photon subtraction is carried out 
using on/off detectors, the resulting state is generally a 
mixed state. Evaluating the entanglement of such a state 
is far from trivial. 

We have applied the negativity and the logarithmic 
negativity, which are known to be entanglement monò- 
tones, to this problem. We have compared the gain in 
terms of these measures with the one in terms of cer- 
tain protocol-specific operational measures, namely the 
teleportation hdelity and the mutual information of the 
dense coding. 

In the asymptotic limit T — > 1, one would have an 
ideal singlo photon subtraction even with on/off detector 
and hcnce have an pure non-Gaussian state, although the 
successful probability of event selection approaches zero. 
This ideal non-Gaussian state is always superior to the 
input squeezed vacuum state for all the measures. For 
T < 1, on the other hand, up to a certain point of A, the 
photon-subtracted squeezed states, both the pure and 
mixed states, are superior to the input squeezed state 
with respect to all the measures. Exceeding this point 
A, the non-Gaussian operation brings no gain, and the 
cntangled squeezed state shows the better performancc. 
When A approaches unity, the effect of the non-Gaussian 
operation for entanglement enhancement gets lost. This 
is because the initial squeezed state approaches the maxi- 
mally entangled EPR state, as A— >1, whose entanglement 
cannot be enhanced by any physical process. 

While the operational measures have clear physical 
meanings, they directly depend on the input state char- 
acteristics, and the evaluation based on them vary for the 
protocols. The (logarithmic) negativity is, on the other 
hand, independent of any such external parameters. This 
quantity reflects the entanglement as an intrinsic prop- 
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erty of the state of interest, namely the final output state. 
We have found that whenever the enhancement is seen 
in terms of the operational measures, the negativity and 
the logarithmic negativity are also enhanced. For the 
dense coding scheme, in particular, the upper limit of A 
below which the non-Gaussian operational gain can be 
seen, A D (pure state case) and Ap (mixed state case) 
approach the ones measured by the logarithmic negativ- 
ity Al N and A^f N , respectively, as the modulation signal 
power f3 gets smaller. It would be interesting to investi- 
gate whether the intersections Ad in the dense coding as 
well as At for the teleportation are identical to Aln or not 
after the operational measures are further optimized with 
respect to all the possibilitics for cxternal paramctcrs of 
input states. In other words, the logarithmic negativity, 
not containing any additional parameter, might give the 
universal upper limit for the region of A, where the gain 
by the non-Gaussian operation can be seen. 

The operational meaning of the logarithmic negativ- 
ity has not completely been clear yet. The relationship 
between the logarithmic negativity and other entangle- 
ment measures, such as the entanglcmcnt of formation, 
has not yet been fully elucidated. Restricted to the case 
of symmetric Gaussian entangled states only, they always 
indicate the same ordering [36( . Generally, however, they 
could give the different ordering for a few cases [3(| |33 • 
Our fmdings will provide some insight into further studies 
of the evaluation of the performance of the non-Gaussian 
operations. 

For practical application of our results to the labora- 
tory experiment, various imperfections should be consid- 
ered, such as the limited quantum cfficiency and nonzero 
dark count rate of the photodetector, and linear loss in 
the optical paths. Thcn the logarithmic negativity eval- 
uation in consideration of these imperfections is the im- 
portant issue. However, those effects cause the complex 
mixing between two modes A and B. Therefore, it is not 
obvious whether the partially-transposed density matrix 
P^q can be split into series of submatrices, as in the case 
of the ideal setup. This would make the analysis much 
harder even by numerical simulation. Such analysis for 
more practical situations is a future problem. 
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APPENDIX A: DERIVATION OF AND itòTI) 

Let us describe the non-Gaussian pure state in terms 
of coherent states, which is suitable for inner products 
with both the Fock state basis and continuous variable 



basis. 



(na 



( x \ a ) = -j]= ex P 
J/n 



x 2 r- a 2 \a\ 2 



(Al) 
.(A2) 



First, we consider the average fidelity of CV telepor- 
tation. With the coherent state basis, the non-Gaussian 
pure state l|l·l·|l is expressed with an auxiliary parameter 
Mi 

I^ng)ab 

VT - A 2 " 



x exp 



d 2 a I d 2 (3RaP 
i±^(H 2 + |/3| 2 ) + Aa*/r 



Ta) A \VT/3) 



rVT^J 2 d 



,(i) d^i 



d 2 a / d 2 f3 



dct 



x exp 



1 + R 



{\a\ 2 + \l3\ 2 ) + \a*l3* + ^ap 



\VTa)A\VT(3} B 



(A3) 



H=0 



where Pjj! is given by (|12|l . and the auxiliary parameter 
\x should be set to zero after all integration and diffcr- 
ential operations. The inner product between the input 
state |ao) and unnormalized output state after telepor- 
tation operation is 



(«o|V p ng P) IVw)b 



(«oIV4g P) ^b(0 TA(H(x,p)\ao) T \ilkt) 



(1)' 



AB 



RVT^J 2 d 1 



2wP 



(i) dfi 1 — \u 



x exp 



dct 



- 1- 



AT 
1 - A^í 



2 XP 



fj,=0 



where 



x + ip 

Q = "o = ç - «o- 

Therefore the (x,p) component of the fidelity is 

NG r NGl X 'W 



(A4) 



(A5) 



iï 2 (l-A 2 



d 2 



i 



27rP dct 9 ^ ld ^ 2 ~ Vl)(l - AM2) 



x exp 



AT 



AT 



1 — A/ii 1 — \U2 



Aíl=P2=0 

(A6) 
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By integrating Eq. (|A6|> with respect to x and p, we above, we obtain the homodyne probability distribution, 
obtain the average fidelity, 



/oc poo 
dx / dpP^ ] F^(x,p) 
-oo J — oc 



(1-A 2 )A 2 T 2 

2^(1 -W) S 



\ 2 T 2 



(A7) 



which is independent of the parameter ao- 

Then, we consider the mutual information of CV dense 
coding channcl. 



AB (U(x,p)\Ü A (x s ,p s )\ V'ng ) ab g — ,í>s {x ~ Xs) 
d 2 a í d 2 f3 



i l-X 2 



exp 



(\a\ 2 + \f3\ 2 ) + \a*/3* +Taf3 



x 2 £' 2 

-Vrç'*a- Vrç'p- — + ^- 



Ra(3 



l-X 



z7r ^det 



exp 



(H 2 + |/?| 2 ) + Aa*/T +//a/3 

-Vrta-Vrç'p-Ç + Ç 



1 1 - X 2 R XT 



XT 



2(1 - XT) 



lí 



/|2 



x exp 



XT 
2(1 - AT) 



ab (ü(x, p) | E/a (a; s , p s ) | ^jjg) al; 



1 - X 2 X 2 T 2 



2nP (r>(l -XT)*\T 



R 



XT 



2(1 - AT) 



li' 



/|2 



x exp 



1 - A 2 A 2 T 2 



1 + AT 
2(1 - AT) 



R 



li 



/|2 



2 7 rFÍ 1 t )(l-AT ) Hr 



AT 



2 £2 



2AT a 



1 + AT/ d/J," 2 í + XTdfi" 



exp 



, 1 + AT 
2(1 - AT) 



li' 



'|2 



(All) 



where the auxiliary parameter fi" should be set to unity 
after all differential operations. With Eq. (|Al·l() . the 
components of channel matrix can be calculated. For 
example, 



p ws (&oo|aoo) 

> p OC 

dx / dpP^\x,p\V2f3,V2P) 
Ja 



o 



1 + erf i 



1 + AT 
1 - AT 



M "/3 



M "=i 



(A12) 



,(A8) where 



where we use the relation 

Ü A (x s ,p s )\Tl(x,p)) AB = e^ p ·< x - x ·>|II(a;-;r, ) p- í >,))AB, 

(A9) 

and 



, (x — + i(p — «„) 
v2 



The auxiliary parameter //(= /i + T) should be set to 
T after all integration and differential operations. From 



2V 



(1 - A 2 )A 2 T 2 



^dct(l — AT) 2 (1 — A 2 T 2 ) \T 



R 



XT 



d 2 



2AT d 



1 + XTJ d/i" 2 1 + XTdfi" 



(A13) 



containing the differential operation with respect to the 
auxiliary parameter /i". Other components can be de- 
rived similarly. With these results, we can obtain the 
mutual information JS7J. 
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APPENDIX B: DERIVATION OF (HI) AND ll62l) 



where 



Similar to the Appendix ^ we describe the non- 
Gaussian mixed state with coherent basis, 



PNG 



1-A 2 
vr 4 Pdct 



exp 



dfai J d z a 2 J <f /?i J d z (i 2 

^(M 2 + M 2 + |/?i| 2 + |/3 2 | 2 ) 



Ta 



Tan 



+A(aí/3* + a 2 (3 2 ) 
1 

T/3i) B (\/T/3 2 



ex J-^^(\ ai \ 2 + \a 2 \ 2 + 1(3^ + 1(3^) 



+A(a*/3* + a 2 (3 2 ) + + JjPiPi 

\VTa 1 )A(VTa 2 \®\VT0 1 )B(VT02\, (Bl) 

where Pdot is give by Eq. ffify. and 71 = R, 70 = 0. 

With the representation i|Bl|) . the unnormalized state 
after teleportation operation is 



PnG ( X >P)$G ( X >P) 

= D B (0TA<n(a;,p)| 
1 



.(in) „ ~(E) 



|n(o;,p))TA^(O ï 



E 

i,i=0 



exp 



1-A 2 
1+ü 



rf 2 /3i / rf 2 /3 2 



-xVtqpi - \VTQ*[3 2 - \Q\ 2 

+^fm-ePi+eí32-m) 

\Vt^+0b(VT(3 2 +^\, (B2) 

where Q and £ are similar to the definition in Eq. I|A5() . 
Therefore, the (a;,p) component of the fidelity is 



^ng (EjP)^ng(íc,p) 
1 



1 



1-A 2 



i,j=0 

x exp 



27rPdet 1 - A 2 7i7j 



-fCi 



(x- Vïa^Y + (p - V2a£ 



(B3) 



l-AT-A 2 77 J --¥ : (7, + 7 3 -) 
= 1 Ï2 

1 (r)^. (i) 

and «o — a o + • 

By integrating Eq. (|B3(1 with respect to x and p, we 
obtain the average fidelity, 



•Fng(A) 



yt-ïf+iJ: 

.4-„ ' 2T dct 1 

h3=0 
1 

E(-ir +j ^(A). 



1-A 2 



-AT-A 2 7l7j - V-(7i+7i) 



(B5) 



For the mutual information of dense coding, we calcu- 
late the homodyne probability distribution, 

Png ( x M x siP») 
= ab(ÏI(x,p)\Üa{x s ,Ps)pngÜI(x s ,Ps)\U(x,p))ab 
1 1-A 2 



27rP det (1 - AT) 2 - A 2 7l7 , 

1-A 2 (T + 7)(T + 7 J ) |CÍ 



x exp 
1 



2{(1-AT) 2 -A 2 7i7i } 



li' 



= J2(-iy + ^(xM^Psh 

i,j=0 

where Ç s and £' are given in Eq. (|A10(i . 



(B6) 



With (|B6fl . the component of channel matrix can be 
calculated. For example, 



^ng ( b oo|ooo) 



OO />00 



dx / dpP^ix^y/ïfryfip) 



E 



(-1) 



1-A 2 



- 47> det l-A 2 (T + 7í )(T + 7,) 



= Eí-^^ill+erf^)] 2 

i,í'=0 



(B7) 



Other components can be derived similarly. With these 
results, we can obtain the mutual information (|62|) . 
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